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There  is  great  interest  in  multiple  input-output  (mio)  feedback  systems, 
for  obvious  reasons.  A  great  deal  of  significant  work  (too  numerous  to  list 
but  Wonham  and  Morse  1972,  MacFarlane  1973,  Wang  and  Davison  1973,  Rosenbrock 
1974,  Porter  and  D'Azzo  1978  are  representative  and  include  bibliographies)  has 
been  done,  primarily  in  the  realization  and  properties  of  the  closed-loop  input- 
output  relations,  under  the  constraint  of  a  feedback  structure  around  the  known, 
fixed  mio  "plant."  There  has  been  notable  work  done  with  uncertain  inputs,  but 
again  only  with  fixed,  known  plants.  Of  course,  plant  uncertainty  is  always 
implicit,  if  only  because  of  the  usual  approximations  required  to  obtain  a  linear 
time-invariant  (fcti)  model. 

In  any  case,  there  does  not  exist  as  yet  any  "quantitative  synthesis" 
technique  for  the  mio  problem  with  significant  plant  uncertainty,  even  for  the 
linear  time-invariant  case.  By  "quantitative  synthesis"  is  meant  that  there 
are  given  quantitative  bounds  on  the  plant  uncertainty,  and  quantitative 
tolerances  on  the  acceptable  closed-loop  system  response.  The  objective  is 
to  find  compensation  functions  which  guarantee  that  the  performance  tolerances 
are  satisfied  over  the  range  of  the  plant  uncertainty.  In  "quantitative 
design,"  one  guarantees  that  the  amount  of  feedback  designed  into  the  system 
is  such  as  to. obtain  the  desired  tolerances,  over  the  given  uncertainty  range. 

In  other  designs,  the  amount  of  feedback  may  be  more  or  less  than  necessary-- 
it  is  a  matter  of  chance.  The  practical  experienced  designer  may  find  the 
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QUANTITATIVE  SYNTHESIS  OF  UNCERTAIN  MULTIPLE 
INPUT-OUTPUT  FEEDBACK  SYSTEM 


Isaac  Horowitz  % 

ABSTRACT  /h  *■  J  ' ~ 

There  is  given  an  n  input,  n  output  plant  with  a  specified  range  of 
parameter  uncertainty  and  specified  tolerances  on  the/ri^  system  response  to 

p 

command  functions  and  the  n  response  to  disturbance  functions.  It  is  shown 
how  Schauder's  fixed  point  theorem  may  be  used  to  generate  a  variety  of 
synthesis  techniques,  for  a  large  class  of  such  plants.  The  design  guarantees 
the  specifications  are  satisfied  over  the  range  of  parameter  uncertainty.  An 
attractive  property  is  that  design  execution  is  that  of  successive  single¬ 
loop  designs,  with  no  interaction  between  them  and  no  iteration  necessary. 
Stability  over  the  range  of  parameter  uncertainty  is  automatically  included. 

By  an  additional  use  of  Schauder's  theorem,  these  same  synthesis 
techniques  can  be  rigorously  used  for  quantitative  design  in  the  same  sense 
as  above,  for  nxn  uncertain  nonlinear  plants,  even  nonlinear  time-varying 
plants,  in  response  to  a  finite  number  of  inputs.  J.  :  v 
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flatter  approach  sufficient.  However,  a  scientific  theory  of  feedback  should 
/  certainly  include  quantitative  design  techniques. 

'In  this  paper  it  is  shown  how  Schauder's  fixed  point  theorem  can  be 
used  to  generate  a  variety  of  precise  quantitative  mio  synthesis  techniques 
suitable  for  various  problem  classes .,^_An^outstanding  feature  of  each  synthesis 
procedure  is  that  it  consists  of  a  success ioh'Tjf  direct  (no  iterations 
necessary)  single-loop  design  steps,  furthermore,  by  a  second  use  of 
Schauder's  theorem,  the  techniques  are  rigorously  applicable  to  quantitative 
synthesis  of  nonlinear  uncertain  mio  feedback  systems.  This  paper  concentrates 
on  existence  proofs  but  a  2  x  2  example  is  included. 


1 . 1  Preliminary  Statement  of  a  Linear  Time  Invariant  MIO  Problem 


In  Fig.  1,  P  =  (p^.(s)]  is  a  n  x  n  matrix  of  the  plant  transfer  functions 
in  the  form  of  rational  functions,  each  with  an  excess  e..  >  0  of  poles  over 

zeros,  and  with  a  bounded  number  of  poles.  The  p.,(s)  are  functions  of  q 

'J 

physical  parameters,  with  m  an  ordered  real  q-tuple  sample  of  their  values. 

M  =  {ml  is  the  class  of  all  possible  parameter  combinations.  The  elements  of 
the  n  x  n  lti  compensation  rational  transfer  function  matrices  F  *  [f^j(s)], 

G  =  (g..(s)j  are  to  be  chosen  practical  (each  with  an  excess  of  poles  (*er 

*  J 

zero).  They  must  ensure  that  in  response  to  command  inputs  the  closed-loop 
transfer  function  matrix  T  =  [tuy(s)]  (of  c  =  Tr)  in  Fig.  1  where  c,  r  are  the 
n  x  1  matrices  (vectors)  of  system  outputs  and  inputs,  respectively,  satisfy 
conditions  of  the  form 


0  <  Auv(“)  i  IV^)!  *  Bu> 
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If  the  t-uv(s)  have  no  poles  or  zeros  in  the  right  half-plane  (are  stable  and 
minimum-phase),  then  tuy(s)  is  completely  determined  by  |tuv(ju))j,  so  (1) 
suffices  (Bode  1945).  It  has  been  shown  (Horowitz  1976)  that  time-domain 
tolerances  of  the  form 

«?(t)  s  <  »J(t) 

dt 

v  =  0,  1 .  n,  any  finite  number,  can  be  satisfied  by  means  of  tolerances 

like  (1)  on  |c(jui)|,  where  c(s)  =  Jc(t).  The  writer  finds  it  much  more 
convenient  to  develop  the  synthesis  theory  in  the  frequency  domain,  and  the 
above  proves  its  sufficiency  for  time-domain  synthesis. 

This  presentation  concentrates  on  the  command  response  problem,  but  the 
same  ideas  can  be  used  to  handle  the  quantitative  disturbance  response  problem 
under  plant  uncertainty,  as  will  be  shown  in  Sec.  6.  The  constraints  on  the 
plant  and  the  specifications  are  introduced  as  needed,  in  order  to  clarify  the 
reasons  for  their  need. 

2.  DERIVATION  OF  SYNTHESIS  TECHNIQUE 
o 

In  Fig.  1,  there  are  available  n  loop  transfer  functions  in  L  * 

2 

(s)J  »  PG,  and  n2  f..  in  F  for  satisfying  the  tolerances  (1)  on  the  n 
tjj.  But  in  the  expansion  of  T  *  (t.j(s)]  =  (I  +  L)  LF,  each  taj)(s»m) 

(mefO  is  a  function  of  all  the  ^(s.m)  each  uncertain,  resulting  in  very 
complicated  expressions  for  tflb  and  making  direct  quantitative  synthesis 
seemingly  impossible— at  least  so  far  unsuccessful.  The  objective  here  is 
to  convert  each  tab(s,m)  design  problem  into  an  equivalent  single-loop  problem 
with  uncertainty.  This  is  done  for  each  tab,  by  lumping  all  the  other  inter- 


acting  t..  variables  into  an  'equivalent  disturbance',  as  follows. 

In  Fig.  1,  c  =  PG(Fr  -  c),  so 

(P'1  +  G)c  =  GFr.  (2) 

Hence,  the  following  restriction  on  P: 

(PI):  A(s)  -  determinant  P(s)  f  O.VmeM. 


Let  r  f  0  and  r.  =  0,  i  j*  v,  so  the  resulting  c  {s)  =  t.  (s)r  .  Let 

VI  J  JV  Y 


p_1  * 

The  uth  element  of  (2)  is  then 


(3) 


rv(s)  tI,  <Pui  *  9uiJtiv  =  ? 

To  simplify  the  presentation,  we  take  gy.  =  0  for  u  f  i  (although  in  practice 
it  may  be  useful  not  to  do  so).  Then  letting  ry(s)  =  1,  the  last  equation  can 


be  written  as 
1 

1 
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(4a) 

(4b) 


This  corresponds  precisely  to  the  single-loop  problem  of  Fig.  2,  with 

d  =  1/P  Of  course,  the  t.  in  d  of  (4b)  are  not  known  but  the  bounds 
*uve  uu  iv  uv 

(1)  on  |tjy|  are  known^ generating  a  set  Dyy  =  (duy>.  We  define  the  extreme  duy 
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Suppose  we  can  find  9uu(s)  and  fuv(s),  such  that  in  the  notation  of  (4,5) 

0  <  Kvl  1  lTduvJlduve>c  [Auv’  Buv]*^meM  •  <6) 

Then  the  magnitude  of  the  right  side  of  (4a)  e[Auv>  Buv]  for  all  mcM  and  for 
all  possible  combinations  of  t^v  (i  f  u)  which  satisfy  (1).  Suppose  this  is  so 

V  u,v  combinations,  and  the  other  Schauder  conditions  of  Sec.  2.1  are 
satisfied.  Then  Schauder's  fixed  point  theorem  can  be  used  to  prove  that 
these  same  n  guu  and  n  f.y  are  a  solution  to  the  synthesis  problem  (1). 

2.1  Application  of  Schauder's  Fixed  Point  Theorem 

This  theorem  states  that  a  continuous  mapping  of  a  convex,  compact  set  of 
a  Banach  space  into  itself,  has  a  fixed  point  (Kantorovich  and  Akilov  1964). 

O 

We  define  the  Banach  space  to  be  the  n  C[0,“)  product  space  denoted  here  by 

C(n  ),  with  norm  =  Z  individual  sup  norms.  C(0,°°]  is  the  Banach  space  of  real 

continuous  functions  f(u>),  we[0H  with  j|f||  =  sup  |f(w)|.  The  convex  compact 

2  ^ 

set  in  each  of  the  n  C[0,<»J  is  taken  as  the  acceptable  set  of  |tuv(jw)|  satis¬ 
fying  (1),  denoted  by  {he(w)}  =  Huy.  Additional  constraints  have  to  be  assigned 
to  the  hfi  (ca)  in  order  that  each  Huy  set  is  compact  and  convex  in  C[0,«°).  These 

constraints  have  been  justified  in  detail  in  (Horowitz  1975)  and  are  therefore 

2 

only  summarized  here.  If  each  set  is  convex  and  compact  in  C[0,~)t  their  n 

2  2 
product  set  denoted  by  H(n  )  is  convex  and  compact  in  C(n  ). 

Constraints  on  Huv  =  (h(a)))  uv 

1.3  continuous  functions  Auy(u),  Buv(w)  with  properties  of  (1)  as 
bounds  on  h(u) 

2.  h'(w)  is  uniformly  bounded:  3  3  |h'(w)|  <K,  V  h.w 


3.  h(w)->o  as  ur*»  in  the  form  k/we,  e  a  fixed  finite  number  >3  to  allow 

at  least  one  excess  of  pole  over  zeros  for  the  elements  of  F,G,P  in  Fig.  1. 
These  constraints  guarantee  (Horowitz  1975)  that  h(u>)  can  be  taken  as  the  magni- 

A 

tude  of  a  function  h(s)„  -  which  has  no  zeros  or  poles  in  the  interior  of 

S~JU) 

A 

the  right  half-plane  or  on  the  jw  axis.  Arg  h(jta)  is  obtained  from  h(m) 

by  anyone  of  a  number  of  Bode  integrals  (Bode  1945). 

o  o 

An  element  of  H(n  )  consists  of  n  positive  functions  on  [0,«),  h.  (to)’. 

.  I  K 

Using  any  appropriate  Bode  integral,  find  the  associated  phase  function  denoted 

A 

here  by  arg[hik(u)] ,  giving  the  minimum-phase  stable  function  h^(s), 

A 

hi|c(j">)  =  hj|c(“)  +  J  atg[h.k(w)].  For  future  use,  denote  this  sequence  of 
operations  whereby  h(u>)  is  transformed  into  h(jw),  as  the  "Bode  transformation" 
B(h(w) ) .  Define  $  on  H(n^)  by 


*  =  **ir*12*  ••**'l'nn}:  H^2^*  *uv^hirh12’  **•»  hnn) 


if  -  y  p  . 

'uu  uv  ^4  ui 


B(hiy(«)) 


P  (1  + 
uu(  Puu' 


(7) 


using  for  P^,  Puu  any  specific  fixed  mcM.  (Note  the  similarity  of  (7)  to  (4a,b)). 

In  Appendix  2,  it  is  shown  that  guu,  fuy  can  be  found  such  that  *  maps 
2 

H(n  )  into  itself.  It  is  also  necessary  to  prove  $  is  continuous,  as  follows. 

0  is  o  continuous  mapping 

2 

is  continuous  if  each  of  its  n  components  is  continuous.  The  first  step 

A 

in  each  mapping  is  B(h.  (w))*1  h^tju).  In  (Horowitz  1975,  Sec.  Ill)  It  is  proven 
that  the  step  h..y(w)+  arg  h^y(u)-0^y(uj)  is  continuous  in  the  C[0,»)  norm.  Hence, 
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the  mappings  h.v(w)->-  h.v(o>)  cose.jv(u))~^v(m),  h^(us)-*  h.v(to)  sin  e^v(w) 

A  O 

“X-V(w)  are  continuous.  The  denominator  of  (7)  is  a  constant  on  H(n  ),  and 

so  are  g  f  and  the  P  .  in  the  numerator.  Thus,  the  numerator  has  the  form 
suu  uv  ui 


Num.  =  |Ka  +  oKfa  -  ?(C.  +  jD.)$.(u)  +  jX.(aj)),  j  =  /-l, 

all  other  terms  real  and  only  thej^,  mappings  on  H(n2).  Infintesimal  changes 
in^v- »  X.  clearly  result  in  similar  change  in  Num.,  so  Num.  is  continuous  on 
H(n2)  and  so  is  each  ij>uv  of  (7)  and  hence  4>.  The  conditions  in  Schauder's  theorem 
are  satisfied,  so  $  has  a  fixed  point. 

This  means  3  a  set  of  h^(u)  denoted  by  h^j(w),  3 


if  .  . 

huv(b,)  * 

' 

9uufuv  Puihiv 

P  (1  +  ?MSL) 
Kuuu  P  ; 
uu 

u»Y  =  1  y  •  •  •  >  J 

n,  where  h*v(j<»)  = 

We  would  now  like  to  deduce  from  (8),  that 


(8) 


B(huyU)) 


huv(ju) 


<Ufuv  -  Ju  PuiC(j“1 

•W  *  iN 

uu 


For,  if  (9)  is  true,  then  by  letting  huv(ju)  =  tuv(ju),  we  have  recovered  (4) 
and  the  n2  h*v(ju>)  are  a  solution  to  the  mio  problem  for  that  specific  meM. 
The  solution  is  unique  if  every  building  block  in  the  mio  system  has  a  unique 
output  for  any  given  input,  which  is  a  very  reasonable  condition.  This  makes 


o 


it  unnecessary  to  prove  that  there  are  no  transitions  from  (8)  to  an 
expression  similar  to  (9)  but  with  right  half  plane  poles  and/or  zeros.  Since 
m  is  any  element  of  M,  this  is  true  for  all  meM  (of  course  with  a  different 


set  of  hyv  for  each  m). 


The  step  from  (8)  to  (9)  is  a  crucial  one  and  must  be  justified  with  great 
care.  Given  an  analytic  function  ${s),  there  is  an  infinitude  of  *(s)  such 
that  |4>(jw)j  =  l^(juj)|,  w  e[0,“],  e.g. 


i»(s)  =  <fr(s) 


(1  -  t i s )  (1  +  t2s) 
TT~+  Tjs)  (1  -  t2s) 


But  4>(s)  /  <l>(s)  even  though  |$(ju)j  =  |o(jw)|*  But  suppose  v/e  know  from  other 
sources  that  4»-| (s )  has  no  right  half  plane  zeros  or  poles,  then  given 

|*1(jt»)|  =  M(u>)  a  magnitude  function  which  is  Bode  transformable,  we  can 
conclude  that  ♦j(jw)  =  B(M(w))  =  M(jw).  Hence,  to  justify  (9)  we  must  prove 
that  the  expression  inside  the  vertical  bars  in  (8)  has  no  right  half-plane 
zeros  or  poles.  The  pole  part  is  easy,  because  1  +  9UU/PU(J  is  obviously 
designed  to  have  no  right  half-plane  zeros;  certainly  guu,  fuv  won’t  be 

A 

assigned  any  such  poles;  v(s )  doesn’t  have  any  by  definition,  and  P  .  is  not 
allowed  any  such  poles--see  Sec.  3.1.  To  prove  the  zero  part,  note  that  from 
(6)  and  Rouche's  theorem,  the  number  of  zeros  of  the  right  side  of  (9)  in  the 
right  half-plane,  equals  such  number  of 


g  f 
Muu  uv 


r  (i  +  ba 

UU'  P 

uu 
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which  is  easily  made  zero  in  the  single-loop  synthesis  steps  (if  PU(J  has  no 
right  half-plane  poles,  a  condition  necessary  for  other  reasons — see  Sec.  3.1). 
Thus,  the  expression  inside  the  bars  in  (8)  has  no  right  half-plane  poles  or 
zeros,  justifying  (9).  This  is  a  very  valuable  result.  The  problem  of 
stabilizing  a  highly  uncertain  n  x  n  mio  system  is  automatically  disposed  of  in 
the  synthesis  procedure,  which  is  furthermore  one  of  designing  n  single-loop 
transmission  functions. 

Tt.  is  worth  noting  that  even! if  the  above  proof  was  not  available,  it 
would  not  be  disastrous  for  this  synthesis  theory.  It  would  only  be  necessary 
to  guarantee  that  at  one  meM,  the  system  is  stable  and  minimum-phase.  For 
then,  this  would  be  soVmeM,  because  by  the  continuity  of  the  poles  (and  zeros) 
with  respect  to  the  parameters,  the  right  side  of  (8)  would  have  to  be  infinite 
(zero)  at  some  w,  in  order  that  for  some  meM  the  system  should  be  unstable 
(or  have  a  right  half-plane  zero).  However,  the  synthesis  procedure  by 
definition  precludes  this.  And  it  is  a  relatively  easy  matter  to  guarantee 
the  desired  conditions  at  one  meM.  . " 

3.  CONSTRAINTS  ON  MIO  PLANT 

The  above  results  hinge  on  our  ability  (a)  to  find  guu  and  f  to  satisfy 

(6)Vu >,  all  u,v  pairs  and  all  me,H  (b)  that  each  equivalent  single- loop  design 

is  stable  and  minimum-phase  V  meft.  These  lead  to  constraints  on  the  mio  plant, 

obtained  by  applying  single-loop  design  theory  to  achieve  (a,b).  Appendix  X 

gives  an  existence  theorem  for  single-loop  design.  The  first  part  of  the  design 

(see  Appendix  A3)  gives  bounds  on  the  nominal  loop  transmission  which  is 

q  /P  of  (4a),  where  P...,-  is  the  ^nominal'  associated  with  a  nominal  mfte|(. 

*uu  uuo  uuo  o 


These  bounds  must  be  satisfied  in  order  that  a  specific  system  transfer 


function  tuy  satisfy  (1).  Here  guu/Puuo  is  used  for  n  tuy  (v  =  1 . n) 

functions.  It  is  proven  in  A3,  that  a  9ul/puu0  can  be  found  which  satisfies 
the  conditions  for  all  n  t  functions. 

For  example,  consider  t^  at  oi  =  and  suppose  ■  -9.  Bul tco, )  » 

1.1  in  (1).  We  could  split  this  range  [.9,  1.1]  into  say  [.95,  1.05]  for 
and  .05  for  T^d^  in  (4),  using  d^  of  (5)  for  d^.  The  technique  in  A3 
or  better  (Horowitz  and  Sidi  1972),  is  then  used  to  find  a  bound  on  9uu(j“])* 
Here,  we  note  a  tough  constraint.  Sooner  or  later  in  uj,  l9uu(jw)|  must  become 
very  small  with  1  +  9uU/Puu  *»■  1  and  then  in  (4a) 


uv 


g  f  -  d 
auu  uv  uv 


uu 


00) 


and  in  (7),  ^uv  -+  the  numerator  of  its  right  side  divided  by  Puy.  How  (4a,  5,  6) 
in  general  require  that 


l^uvLax  >  ^lTduv^uve^ 

But  l^jJmax  =  Buv  and  at  high  ^re9uencies 


01) 


rLlhnl|B.vl 
lTduvduvel  *  [E^l  ‘ 

To  see  what  this  leads  to  take,  for  example,  n  =  2  so  that  the  above  applied 
to  v  =  1 ,  u  -  1 ,2  gives 


> 


2 1^12 !B21 


J21 


yP21iBn 

1^22  * 


11 


requiring 


1  HV21I 
lPllP22l 


as  to  -y  to 


02) 


Thus,  a  constraint  on  P  is 


^  0)^,  ^  for  w  >  l^nP22 1  >  ^P12P21 


(13) 


It  is  known  that  as  s  +  ®, 


ij 


hL 

•u  ’ 


so  the  above  becomes 


1^11^22 1  4 1^12^21 1 


ell+e22 
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e12+e21 
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If  the  uncertainties  in  the  k^  are  independent  and  e^  +  =  e^  +  ^2^  > 

this  becomes 


*cllmin*c22min  >  ^^max^lmax' 


04) 


There  is  an  important  problem  class  for  which  the  inequality  is  less 
harsh.  This  is  the  "basically  noninteracting"  class,  where  one  ideally  desires 
t. .  =  0  for  i  i  j,  but  because  of  uncertainty  accepts  A. .  =  0,  |t.-|  <  B. . 

IJ  l  J  I J  I J 

for  i  f  j,  in  (1).  Also,  one  doesn't  care  if  t - . ( i  f  j)  is  nonminimum-phase. 

*  J 

Condition  (6)  then  applies  only  to  u  =  v.  The  f  (u  t  v)  are  set  equal  to 
zero  and  (13)  becomes 

3  lPllP22 1  >  ^P1?P?1  I  *ne^*  u  >  wh’ 


05) 


It  is  desirable  to  ease  inequality  (13)  in  the  general  case.  Note  that  (6) 
can  be  satisfied  over  any  fini te  w  range  by  making  |1  +  g  /p  |  large  enough. 

Thus,  as  previously  indicated,  one  can  split  the  [Auv»Buy]  tolerance  so  that 
1TUV!  >  lTduv^duve!’^meMj  e‘g-  assi9n  lTuvl  £  [E  -  e,  E  +  e]  with 

^  =  ^uv  +  BUV)/2’  2t  <  Buv  '  \v  and  the  balance  (Buv  "  -  2e)/2  is  assigned  to 

Tduvduv  ^a)*  Bu^  |1  +  9UU/5UU |  must  then  be  made  large  enough  to  satisfy 
uhe  resulting  requirements,  and  it  can  for  any  finite  w  range.  The  trouble  is 
that  9uu  must  be  allowed  to  ->■  zero  as  u  +  leading  to  (13),  etc.,  if  we 
insist  on  (6).  We  could  ignore  (6)  at  large  u,  say  for  to  >  Wu,  with  Wu  as 

n  H 

large  as  desired  but  finite,  letting  j-ruv|  «  |Tduv(|duve|  for  u  >  ^  Then 
for  u  >  w^,  (11)  is  replaced  by  the  weaker 

*  lTduvduvel  (16) 

and  for  n  =  2,  (13)  is  then  replaced  by 

(P2b)'-  3  «v  ^  for  u  >  °V  ^P1  lP22 ^  >  IP12P21 07a) 

-An  important  question  is  whether  (17a)  is  an  inherent  basic  constraint  in  the 
presence  of  uncertainty,  no  matter  what  design  technique  is  used,  or  is  due 
only  to  this  specific  design  technique.  The  methods  suggested  in  (Rosenbrock 
1974,  Owens  1978)  to  achieve  diagonal  dominance,  may  be  helpful  in  satisfying 
(17a),  but  they  would  have  to  be  extended  to  uncertain  plants.  Note  that  in 
Rosenbrock  1974,  Owens  1978),  diagonal  dominance  is  desired  V  w  e  10»°°)» 
whereas  in  (P2b)  ^  required  only  for  u  >  u^. 

For  the  analog  of  (17a)  at  n  =  3,  it  is  found  that  diagonal  row  dominance 
of  P for  0)  >  is  a  sufficient  condition.  The  necessary  condition  can  be 


written  as 


3 


3  V  ?for">“H  lpiipjjl  5  and 

IPnp33l  >  ^  lP12P23^  +  lP13P22^^P22P31 1  +  lP21P32^  07b) 

which  can  be  written  as, 

lPllP22P33l  >  lP11P23P32l  +  ^P12P21P33^  +  lP12P23P3ll 

+  lP13P22P3ll  +  IP13P21P32  |  for  u  >  ujh  .  (17c) 

The  latter  has  the  following  intepretation.  Array  the  matrix  P~^  in  the  usual 
manner,  but  twice -one  under  the  other  as  in  Fig.  3a.  Then  the  terms  on  the 
right  side  of  (17c)  consist  of  the  products  of  the  entries  crossed  by  the 
dashed  lines. 

However,  if  is  so  used,  it  is  no  longer  possible  to  use  Rouche's 
theorem  and  thereby  prove  each  t . .  is  minimum-phase.  But  we  can  still  design 

V  J 

so  that  the  nominal  t. .  are  minimum-phase  and  we  know  from  (6)  that  t.-(ju)  ?  0 

•  J  1  J 

for  o)e[0,uh].  Therefore,  from  the  continuity  of  the  zeros  of  t..  with  respect 
.to  the  parameters  of  the  system,  if  t. .  has  any  right  half-plane  zeros,  they 

*  j 

must  enter  the  right  half-plane  as  shown  in  Fig.  3b,  It  is  unlikely  that  such 
a  zero  which  must  migrate  all  the  way  up  to  jui^,  should  move  back  into  the 
*  significant  control  bimdwith  region  A.  The  point  is  that  if  right  half-plane 
zeros  are  "far-off",  they  have  little  effect  and  the  system  is  "dominantly", 
minimum-phase. 


Rouche  s  theorem  can  still  be  used  if  we  can  guarantee  that  (6)  is 
satisfied  for  a  semicircle  consisting  of  the  segment  and  the  right 

half-plane  half-ci rcumference  of  the  circle  of  radius  Uu,  centered  at  the 
origin.  Then,  there  are  definitely  no  right  half-plane  zeros  of  t. .  in  this 
half-circle,  and  the  system  ’S  "dominantly"  minimum-phase  This  is  quite 
practical  in  the  design  technique  of  (Horowitz  and  Sidi  1972),  discussed  in  A3. 


3.1  Modification  of  mapping  $ 

Note  that  for  the  "dominantly  mini  mum- phase"  and  the  "basically  noninter¬ 
acting"  cases,  the  application  of  Schauder's  theorem  in  (2.1),  Eqs.  (7-9),  etc., 
needs  modification,  because  nonminimum-phase  tuv(ju)  cannot  be  uniquely 
derived  from  |t  (jw)|.  Redefine  h  e  Huy  of  2.1  to  consist  of  an  ordered 

A 

pair:  h(w)  as  before  and  q(w),  the  imaginary  part  of  huy(jw)  with 

A 

h  *  |huv(ju))|;  h  e  Huy  the  same  as  before  but  q(cu)  e  C  [0,«)  with 

0  £  | q (uj) |  <  h(o>).  Constraints  2,3  in  2.1  on  h(w)  also  apply  to  q(ia).  Let 

(HQ)UVCC2  10.-)  denote  the  set  {(h(w),  q(u>))}  with  ||(h,q)||  =  ||h||  +  ||q||. 

Obviously,  (NQ)uV  is  compact  and  convex  in  C  (0,®).  The  extension  to  the 
2 

n  product  set  is  straightforward. 

The  mappings  <|juy  in  (7)  are  redefined.  Each  i|>uy  is  a  pair  of  mappings, 
one  the  absolute  value  as  before,  the  second  the  imaginary  part  with  the 
absolute  bars  on  the  right  removed.  On  the  right  side  of  (7),  Bfh^f&j))  is 
replaced  by  r.v(u>)  +  jqiv(u>),  with  h?y  =  r?y  +  q?y,  (h.y,  qiy)  e  (HQ)iy. 

It  is  necessary  to  prove  that  $  maps  each  element  of  (HQ)uy  into  itself. 


15 


The  proof  follows  immediately  from  that  for  the  minimum-phase  case  --  this 
is  obvious  from  (6),  the  definition  of  duvg  in  (5),  and  Appendices  1,2.  The 
proof  that  $  is  continuous  is  straightforward.  Accordingly,  the  Schauder 
conditions  are  satisfied  and  there  exists  a  fixed  point  which  satisfies  the 
specifications.  Such  specifications,  by  themselves,  would  not  be  good  ones 
because  they  permit  highly  nonminimum-phase  tuv(s).  However,  they  are 
satisfactory  if  it  is  known  from  other  sources  that  t  is  "dominantly 
mini mum- phase". 

3.2  Additional  Constraints  on  P 

Constraints  Al(l)-(3)  in  the  Appendix,  must  be  applied  to  the  1 1 puu» 

since  in  Fig.  (2)  p  =  =  p  of  Appendix.  Al.l  requires  that  there  be 

uve  *uu  ,  . 

no  change  in  the  excess  of  poles  over  zeros  of  p —  =  -r —  where  A  =  det.  P 

ruu  fluu 

and  Auu  its  uuth  minor,  as  m  ranges  over  Also,  that  for  at  least  one 
mcjl,  denoted  by  mUQ,  Puu  has  all  its  poles  and  zeros  in  the  interior  of  the 
left  half-plane.  The  mUQ  can  be  different  for  each  u. 

A1.2  requires  that  1/P  is  minimum-phase  V  me|(,  and  its  zeros  do  not 
get  arbitrarily  close  to  the  jw  axis.  Since  1/PUU  -  A/Auu,  this  means  A  must 
have  no  right  half-plane  zeros.  Hence  the  P^  in  general  have  no  right 
half-plane  poles.  (For  those  who  wish  it,  P  is  restricted  to  be  controllable 
and  observable  V  mefl,  but  these  concepts  are  unnecessary  if  P  is  properly 
formulated  in  terms  of  physical  uncertain  parameters  (Horowitz  and  Shaked 
1975)).  Since  the  p4i  in  P  =  [p..]  are  finite  rational  functions,  the  latter 
part  of  A1.2  is  automatically  satisfied. 


A1.3  for  n  =  2  is  the  same  as  (17),  which  shows  that  (17)  is  a 
fundamental  condition  for  linear  time-invariant  design,  not  an  "extra" 
condition  due  to  our  design  technique,  at  least  for  n  =  2.  However,  (13) 
is  an  "extra"  condition.  Note,  the  extension  of  single-loop  design  to 
disappearing  poles  and  zeros  in  A6  may  perhaps  permit  disappearing  poles 
and  zeros  in  the  mio  plant  functions. 


4.  OTHER  DESIGN  EQUATIONS 


The  previous  design  equations  constitute  only  one  of  many 
design  techniques  derivable  from  Schauder's  fixed  point  theorem.  Only 


two  more  will  be  briefly  mentioned  here. 

Both  are  based  on  the  use  of  a  nominal  diagonal  loop  trans¬ 
mission  matrix.  The  design  obligations  on  the  loop  transmission 
elements  are  then  independent  of  the  way  the  plant  input  and  output 
terminals  are  numbered.  If  G  is  m<ide  diagonal,  such  numbering  is 


important  and  after  one  arbitrarily  numbers  the  plant  input  terminals, 
he  should  try  to  number  the  outputs  such  that  the  main  effect  of  in¬ 
put  i  is  on  output  i.  Manipulation  of  (2)  somewhat  differently  from 
>  Sec.  2,  gives 


kn  i  +  *n/6n 


f„,*„o/<5 


l  v_.t.,/6,0,  etc. 
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=  '2122'  22  i*; 


'2i viV  22 


1  +  *22/6Z2 
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where  V  =  (v.  .]  =  I-^(P)  ,  P  is  the  'nominal' plant  matrix  and  there- 

1 J  U 

fore  fixed,  P  is  the  general  uncertain  plant  matrix,  6..  =  1  -  v... 

r  n  ii 

The  .  are  the  nominal  elements  of  the  loop  transmission  matrix  L. 
Eqs.  (18)  lend  themselves  to  single-loop  design  and  use  of  Schaiider's 
theorem,  pricely  as  did  (4). 

Another  interesting  set  of  design  equations  is  obtained  by  de¬ 
signing  to  control  the  changes  in  t..,  rather  than  t. .  directly.  Let 

1 J  1  J 

Tq  =  lt^ol  be  the  'nominal'  system  transfer  matrix  and  T  =  [t„)  the 

actual  which  is  uncertain,  AT  =IAt..]  =  T  -  T  .  Then  it  can  be  shown 

ij  o 


AT  =  (I+Lf 1  VT,  V  = 


where  are  likewise  the  'nominal'  and  uncertain  plant  transfer 

matrices,  and  L  -  PQ6  =  is  the  nominal  loop  transmission  matrix. 

If  L  is  taken  diagonal,  the  result  is  (n  a  2  for  simplicity) 


*11*11  *  *12*21 
1+4,, 


Vllt12  *  *12*22 
1  + 


and  similar  obvious  ones  for  ,  At^g- 

The  design  problem  is  now  completely  one  of  disturbance 
attenuation,  with  the  disturbances  d^  =  v^tjj  +  v^t^,  etc’»  whose 
range  is  known.  Schauder's  theorem  is  applicable  in  the  same  manner 
as  before.  Note  that  V  represents  the  'normalized'  plant  variation 
matrix.  Eqs.  (20)  appear  to  be  much  simpler  to  use  for  design  (once 


the  At.,  tolerances  are  formulated)  than  (4),  and  their  use  needs  to 
^  J 

be  intensively  researched.  However,  both  for  (18)  and  (20)  the  con¬ 
straints  considered  in  3.,  leading  to  (11-15)  must  be  found,  and 
these  may  possibly  be  tougher  than  before.  Also,  both  a  nominal  2 
and  T  must  be  chosen^which  is  not  goodj because  the  optimum  pairing 
is  not  apriori  known.  However,  the  analogs  of  (14,17)  may  be  more  lenient. 

4.1  Bandwidth  Minimization 

An  important  criterion  for  comparison  of  design  techniques  is 
their  "cost  of  feedback,"  which  we  take  as  the  bandwidths  of  the 
loop  transmission  functions--because  they  determine  the  system 
sensitivity  to  sensor  noise.  Obviously,  quantitative  synthesis 

techniques  must  first  be  invented  before  one  can  turn  to  their  op¬ 
timization  (for  without  such  quantitative  techniques  comparison  is 
possible  at  best,  by  analysis  after  a  specific  numerical  design  has 
been  made).  This  approach  via  Schauder's  theorem  promises  to  generate 
a  variety  of  such  techniques,  and  the  next  step  will  be  optimization. 

5.  DESIGN  EXAMPLE 

The  2  x  2  plant  elements  are  p..  =  k../(l+sA..)  with  correlated 

1 J  I J  I J 

uncertainties,  giving  a  total  of  9  parameter  sets  in  Table  1.  The 
design  was  performed  to  handle  the  convex  combination  generated  by 
these  9  sets  (Figure  6). 


TABLE 


A  "basically  noninteracting"  system  is  desired,  with  the  off-diagonal 
transmissions  specified  in  the  w-domain  |tj2(jw)|,  |t21(ju»)|<0.1  Vw.  The 
diagonal  t^,  t22  bounds  are  identical  and  were  originally  in  the  time- 
domain  in  the  form  of  tolerances  on  the  unit  step  response  shown  in 
Fig.  4a,  b  (which  also  shows  the  design  results  for  those  of  the  9  cases 
which  were  reasonably  distinguishable).  These  time-domain  bounds 
were  translated  into  the  "equivalent"  bounds  on  |t..(ju))|  shown  in 

Fig.  5  (Horowitz  and  Sidi  1972,  Krishnan  and  Cruickshank  1977). 

» 

Familiarity  with  quantitative  single-loop  design  is  assumed 
here  •  One  can  do  a  problem  of  this  complexity  by  hand.  The 
sets  (p..  (ju>)),  called  the  plant  templates,  are  obtained  on  the 

_i 

-1  L  x^'  A 

Nichols  chart.  Some  of  these  templates  of  *  — 


are  shown  in  Fig.  6  at  various  w  values.  The  larger  the  template, 

the  greater  uncertainty  at  that  ui  value.  The  tolerances  on  t  of 

uu 

(4a)  and  Fig.  5  were  divided  between  and  Tdui(d  as  discussed  in 

Sec.  2.  Each  of  these,  in  conjunction  with  the  templates,  leads  to 

g 

bounds  on  the  nominal  loop  transmission  l  =  Some  of  these 

UCOO  £ 

uuo 

boulds  on  due  to  are  shown  as  solid  lines  in  Fig.  7,  i.e., 

it  is  necessary  for  Jlllo  to  lie  above  the  indicated  boundary.  The 
tolerances  on  Tduuduu  lead  to  the  dashed  line  bounds  on  «11q.  No 
attempt  was  made  to  optimize  the  division  of  the  tolerances  between 
Tj -J  and  t^d^.  The  composite  bound  on  must  satisfy  both. 

The  *llo(jw)  chosen  is  also  shown  in  Fig.  7.  There  was  no  attempt 
made  to  optimize  the  i.^;  the  design  was  made  by  hand  quickly,  so 
the  are  larger  than  need  be,  with  the  tolerances  therefore 

satisfied  better  than  necessary— as  seen  in  Figs.  4a,  b.  Optimal 
tii0(ju)  would  lie  on  their  boundariesat  each  u,  so  in  this 
example  there  is  considerable  overdesign. 

Here  we  took 
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v/i  th 


Ao  .75  (1+3. 66s) 
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f,lo  ~  (l+3s)(l+2s) 

The  requirements  on  f]r  (f]2  =  f£1  =  g12  =  921  =  0  here)  were 

found  using  single-loop  design  technique  [15]  as  briefly  explained 

here  in  A4,  and 

f  =  — - —  f  =  — y 

11  1  +  .5s  *  t22  1  +  .33s 

were  found  satisfactory.  The  system  was  simulated  on  the  digital  com¬ 
puter  with  the  results  shown  in  Figs.  4a,  b.  The  t12>  \  tolerances 
were  easily  satisfied  by  the  design. 

While  this  is  not  a  very  challenging  example  of  the  design  techni¬ 
que,  nevertheless  the  uncertainty  is  very  large  and  one  should  consider 
how  quick,  simple  and  straightforward  was  the  design  procedure,  and 
also  consider  what  alternatives  are  offered  in  the  mio  literature. 


There  are  no  other  techniques  available  for  systematic  design  to 
specifications  in  the  presence  of  significant  uncertainty,  which 
guarantee  design  convergence  and  attainment  of  design  tolerances. 


Whatever  present  popular  technique  is  used,  it  would  be  necessary 
to  cut  and  try  and  endeavor  to  understand  the  relations  between  the 
cutting  and  the  results  as  one  continued  to  cut  and  try,  because 
these  techniques  have  no  provision  for  significant  uncertainty.  In 
the  above  design,  one  sweep  was  known  to  be  sufficient  because  the 
plant  and  the  design  tolerances  (w-domain)  satisfied  constraints, 

PI  etc. 

5.  EXTENSION  TO  NONLINEAR  UNCERTAIN  MIO  PLANTS 

Once  there  is  a  quantitative  design  technique  for  linear  time 
invariant  mio  uncertain  plants,  it  appears  at  least  conceptually 
possible  to  extend  it  to  a  significant  class  of  nonlinear,  even 
nonlinear  time-varying,  uncertain  mio  plants.  The  procedure  is  a 
generalization  of  that  used  (based  also  on  Schauder's  theorem) 
in  (Horowitz  1976)  for  single  loop  uncertain  nonlinear  systems. 

The  key  feature  is  the  replacement  of  the  nonlinear  plant  matrix 
set  (a  set  because  of  the  uncertainty),  by  a  linear  time  invariant 
plant  set  which  is  precisely  equivalent  to  the  original  nonlinear 
set,  with  respect  to  the  acceptable  system  output  set.  The  pro¬ 
cedure  is  briefly  presented  for  the  case  where  one  wants  the  system 

with  nonlinear  uncertain  plant  to  behave  like  a  linear  time-invariant 
system  for  a  specified  class  of  command  input  sets. 


It  is  essential  that  the  command  input  sets  represent  a  good 
sampling  of  how  the  system  will  actually  be  used.  For  example,  suppose 
n  =  3  and  in  actual  use  r^,  always  exist  simultaneously  (with 
=  0),  and  appears  by  itself  (with  r]  =  r2  “  0).  Say  there  are 
ten  typical  r^(t)  inputs  and  for  each  typical  r^(t)  there  are  five 
typical  r^(t).  This  makes  a  subtotal  of  50  input  sets,  to  which  is 
added  the  number  of  typical  r3(t)  say  10,  giving  a  class  R  =  {r}  of 
60  sets,  of  which  50  have  the  form  r  =  (r^ ,  r^%  0)  and  r  =  (0,  0,  r3) 
for  the  balance.  Choose  r^eR.  The  family  of  acceptable  outputs  for 
this  input,  is  known  from  the  tolerances  on  t..,  giving  for  that  one 
input  vector  a  family  -H  -  (h>,  h  -  (b^,  h^,  h^).  The  mio  plant  is  re¬ 
presented  by  a  family  (because  of  parameter  uncertainty)  W  of  nonlinear 
differential  mappings 

W  -  {W>,  W  =  (w1,w2,w3)^c1  =  w1(x2,x2,x3>m),  .  .  .  ,  c3  = 
w3(x1,x2,x3,m),  where  the  are  the  plant  inputs,  c.  the  plant  outputs, 
and  m  is  the  plant  parameter  vector  meM. 

Take  a  sample  acceptable  outputtriple  h  =  (h^,  h2>  h3)  and  find 
the  correspond!' ng  plant  inputs  at  some  specific  m&U  (or  in  other  words, 
pick  a  WeW)  and  let  c.  =  h.  and  solve  the  nonlinear  equations  backwards, 
giving  the  input  set  (x1 ,  x2>  x3).  Take  the  Laplace  transforms  ^(s) 
of  X.,  fi.(s)  of  h.  giving  the  vectors  x[s]  =  (x^s),  x2(s),  x3(s)) 

B [si  =  (fi^s),  ....  fi  (s)).  Repeat  for  other  h  samples  in  the 
acceptable  output  set  H,  giving  two  paired  families  of  £[s],  n.[s). 

w 


Select  any  combination  of  three  x[s),  forming  a  3  x  3  matrix  X  and 

A 

corresponding  paired  combination  of  three  n[s],  forming  the  matrix  H 
Set  H  =  PX  and  solve  for  P  =  H(X)  \  P  is  the  linear-time-invariant 
equivalent  of  the  specific  WeW  picked*  with  respect  to  the  specific 
trio  of  acceptable  output  vectors  picked.  Repeat  over  different 
trios.  Repeat  the  entire  operation  over  different  wet1,  giving  a 
class  Pj  =' {P}  which  is  the  linear-time-invariant  equivalent  of  the 
W  family,  with  respect  to  the  class  of  acceptable  outputs  H  for  in¬ 


put  vector  r.p  Repeat  the  entire  operation  for  r^,  .  .  .  ,  r^^  giving 
{?i  }  =  ptota-|  which  is  the  linear  time  equivalent  for  the  nonlinear 
W,  with  respect  to  the  tribe  of  60  families  of  acceptable  output  sets. 


The  equivalence  is  exact  if  the  conditions  for  application  of 


Schauder's  theorem  are  satisfied.  We  now  have  a  linear  time-invariant 


uncertain  mio  problem,  which  let  us  presume  we  can  solve.  If  and  only 
if  we  can  guarantee  the  solution  of  the  latter,  then  the  same  compen¬ 
sation  functions  will  work  for  the  original  nonlinear  uncertain  mio 
plant.  Hence  the  importance  of  quantitative  linear  time  invariant 
design  techniques  (over  and  above  their  intrinsic  importance) — for  they 
enable  the  precise  solution  of  nonlinear  uncertainty  problems. 


The  design  effort  in  the  above  appears  to  be  enormous  but  it  is 
conceptually  straightforward  and  easy.  An  ordinary  control  engineer 
can  implement  it  and  the  digital  computer  is,  of  course,  an  essential 
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tool.  Conceptually  too,  it  appears  possible  to  extend  the  method  to 
obtain  nonlinear  relations  between  inputs  and  outputs  within  specified 
bounds,  despite  large  plant  uncertainty,  even  nonlinear  time-varying, 
as  can  be  done  for  the  single  input-output  case.  The  prospect 

is  fascinating.  Imagine  being  able  to  work  with  the  actual  nonlinear 
equations  of  a  jet  engine,  or  a  chemical  process,  etc.,  include  un¬ 
certainties  in  the  modelling,  even  uncertainty  in  system  order  (see 
Appendix),  and  designing  to  achieve  outputs  within  specified  tolerances 
over  the  given  range  of  uncertainty. 


6.  DISTURBANCE  ATTENUATION 


Let  x  in  Fig.  1  be  a  nxl  matrix  of  distrubances.  The  resulting 
system  output  (with  r  «  0)  is  c  =  (I  +  PG)’1  Px  4  Zx.  Z  =  [t^l.  the  o  x  n 
disturbance  response  matrix.  Bounds  on  Z  are  given  in  the  form 


l2uv(Jto)l  <  buv(u)  *  ^ 

Rewrite  c  =  Zx  in  the  form  (P'1  +G)c  =  x.  Let  xi  t  0  only  for  i  =  v,  so 


(21) 


ci  =  zivxv*  and 


l  <pui  +  9gi)2iv  '  6»  ’  (l'.  u  ■  I1  ’ 


i=l 


*Puu  +  9uu)zuv  =  6v  '  ^  (pui  +  9ui^ziv 

i^u 

i  (pui+9ui>ziv 

_  =  JjLu_ _  (22) 

uv  q 


Let 


‘„ve<“>  *  S“P 
m 


l 

i^u 


u  i 


'ui 


uu 


|biv(uj) 


(23) 


The  gui(w)  (i7u)  can  be  chosen  to  minimize  xuve(o>),  but  for  simplicity  we  shall 
assume  them  zero.  From  (22,23) 


'uv 


M|  <  | 


5Xu  1  *uve, 

0  *  jr^) 

uu 


(24) 


If  1/PUU  satisfies  the  constraints  listed,  then  it  is  obviously  possible 
to  guarantee  lz  (w)|  <  any  finite  number,  no  matter  how  small,  at  any  finite 
w.  Also  it  can  be  made  zero  at  a  finite  number  of  w  values  by  assigning  poles 
to  guu  at  these  values.  Assume  that  guu  can  be  chosen  to  satisfy  (21) 

V  (j  e  [0,°°).  Then  one  can  set  up  the  conditions  for  Schauder's  theorem, 
precisely  as  was  done  in  2.1.  The  set  buv(ui)  must  have  been  formulated  such 

O  O  A 

that  B(n  ),  the  n  product  set  of  the  buV(oj),  is  compact  convex  in  C(nc)j 
o 

analogous  to  H(n  )  in  2.1.  The  analog  of  4>  in  (7)  must  be  formulated 

with  the  modification  of  Sec.  3.1,  inasmuch  as  we  do  not  care  if  the  zuv(s)  are 

nonminimum-phase. 

Conditions  analogous  to  (12-17)  for  n  =  2,  are  obtained  as  follows. 

As  w  ■+  °°,  9UU/PUU  0  so  in  (24),  the  right  side  -*■  its  numerator.  But 
I zuv( jto)  I  £  buv(u»)  of  (21).  Let  u  =  1,  v  =  2  and  then  u  =  2,  v  =  1  and 
obtain  the  necessary  condition  (for  g^  =  =  0), 


As  u  -*■  »  , 


^12^21  <  P 11^22  ’  ^  me:M. 


(25) 


similar  to  (17)  but  here  only  at  because  there  is  no  concern  with  the 
minimum-phase  property.  Setting  u  =  v  =  1,  and  then  u  =  v  =  2  in  (24),  we 
get  the  conditions 


As  to  •»  «°,  bn  >  IpJ— |  =  |pn 


11 


P12P21 

p22 


I  *  ^22  >  I  p22  ~ 


P12P21 

P11 


(26) 


But  in  reality  as  to  ->  «,  c  ->  Px  so  Z  -*  P  and  ->  p^,  z22  -+•  p22.  Hence, 
assignment  of  b^.  (as  to  -><*>)  to  satisfy  (25)  is  no  obstacle,  because  the  buv(to) 
are  upper  bounds  on  the  J  zuv ( joj)  J . 


APPENDIX  1 


EXISTENCE  THEOREM  FOR  SINGLE-LOOP  DESI G N 


The  plant  transfer  function  p(s)  is  uncertain,  belonging  to  a 
set  P  =  {p (s ) }  and  is  imbedded  in  a  two-degree-of-freedom  single¬ 
loop  feedback  structure,  as  in  Fig.  2  (p  in  place  of  p  ).  The 

'uve 

rational  functions  f(s),  g(s)  (replacing  fuv>  g^u  in  Fig.  2)  are 
to  be  chosen  to  satisfy  specified  tolerances  on  the  command  fre¬ 
quency-response  t(jio)  =  and  disturbance  frequency  response 
td(jw)  =  c(jw)/d(jw),  (r,  d,  c  replacing  r  ,  -d  ,  cu  in  Fig.  2). 


A1 .  Constraints  on  P 

1.  p(s)  is  a  rational  function  with  a  fixed  excess  e>l  of 

poles  over  zeros  (this  is  relaxed  later  in  A6,  7).  3  at  least  one 

pcP  one  of  which  is  designated  as  p  ,  all  of  whose  poles  and  zeros 
are  in  the  interior  of  the  left  half-plane. 

2.  At  each  we[0,°°),  3  inf.  |p(ju)|  =  b(w)>0.  3  inf  b(u) 

P  I 

=  bj>0  for  any  finite  interval  I  =  10, ul.  Also,  |pQ|  of  Al(l)  has 
a  sup  on  each  finite  interval  I  =  10, ul,  sup  |pQ|  =  xqj. 

3.  As  s-xo,  p(s)-*-kp/se,  k  e[kj,  kgl  with  «>  k^k^O,  uniformly 
on  p  in  the  following  sense:  For  any  e>0,  no  matter  how  small,  3  ue 
(independent  of  p(s)),  such  that  for  each  peP  there  is  associated  a 


k  efkj,  kg]  so  that 


|ani  <e  and  Arg  |  p(ju)+e^-|<e,  for  u  y  Ufa  • 
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Note  that  Al(l)  permits  changes  in  plant  order,  e.g.,  - 

v/ith  say  a  e  [2,  5],  [0,  3].  Al(2)  dictates  minimum-phase  p(s) 

and  that  the  jw  axis  is  not  a  limit  of  any  sequence  of  p(s)  zeros. 

A1 (3)  requires  a  uniform  bound  on  the  poles  and  zeros  of  all  peP. 

A2.  Tolerances  on  1 1 (.jco)  1  and  |t^(.itj)| 

(1)  0<A(w)  t(jw)  |  ^B(oi)  with  A,  B  e  C  [0,  »),  »b(w)>1. 

A(u) 

3  inf  3(w)  =  3.>1  on  any  finite  I  =  [0,uj]. 

1  Xk 

(2)  3  X  >  1.05,  w,  3  for  ,  ^4  =  This  means  that 

A  A 

in  the  high  w  range,  the  feedback  is  allowed  to  increase  the  sensitivity 
S  =  ^p(jyj/p(ju)’  rat^er  ^crease  it.  In  fact,  as  noted  by 

CO 

Bode,  j  ln|S|da)  =  0  in  any  practical  system,  so  the  decrease 

0 

in  S ( | S | <1 )  achieved  in  the  control  bandwidth  range,  must  be  balanced 
by  | S | >1  in  another  range.  X  can  be  a  large  number,  because  as  or*®, 
e.g.,  suppose  k2/k]  =  10,  who  cares  if  |t(jw)|  ellO"11,  10”7] 

(X  =  10  )  at  very  large  u>. 

(3)  The  tolerances  on  t^jw)  are  in  the  form  | trf( jto) |^Q(o>)>0. 

For  any  I  *  [0,  m],  3  inf  Q(w)  =  3^.  Since  t^  =  P(1  +  pg)"^  = 
pS  of  A2(2),  { Q/p | >1  at  high  frequencies,  so  3  3  for  uxu^, 

Q(<*>)  =  B-j(w)  |p(jo)| ,  B-j >1  - 05. 


stricted  to  the  interior  of  the  left  half-plane  VpeP,  and  minimum- 
phase. 


A3 .  Choice  of  Compensation  Functions 

Let  po(s)  of  A1 (1 )  be  the  'nominal'  plant  with  kQe[kj,  k2l  its 

associated  k  value  of  Al(3). 

P 

Let  £  =  .01  kQ  in  Al{3),  wt  =  largest  of  {w£»  Uy  <odh  It  =  [0,wt] 
kT 


^  =  sup  W, 


•t  \ 


of  A1 (2) .  In  Fig.  2, 


fi. 


t(s)  =  ,  l  -  g  p  . 

1+gp  p  o  3  Ko 

—  +$. 

P  o 


(Al) 


Slip  |t(jaj)| 

“e  want  TnfTtU^T  *  AM  of  ft2(1)-  ™5  is  ach1eved 


in 


P  .  i 

sup 

0+0 

p 

P 

inf 

P 

prt  +  l 

0  0 

p 

<8  of  A2(l ). 
t 


Since 


U0l+Yt 

<  Yt  in  It,  it  suffices  for  It>  if  l*0l>Ytand  JjfJZ  y~  <  P\  » 

0  t  t 


giving  the  sufficient  condition 
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To  satisfy  A2{3)  in  1^,  it  is  necessary  that  |t  |  = 


1+gpj 


P  o 


<0(0)),  which  is  certainly  achieved  if 


l*ol>  sup1p0I  sup|po|  x 
flllt  +inflpl"  = 


1  .  JL\  A 

OIt  IV  bi  '  ■ 1  • 

c  t  I 


od1 


Therefore,  choose 

|*o(ju)i>  larger  of  (|£ot|,  l£od! )=l Aoxl »  In  1 

Next,  we  find  a  bound  on  £q  in  It  =  to  satisfy  A2  in  It. 

From  A1  (3) ,  in  I  =  [u.,  «),  lies  in  the  narrow  sliver  V  in 

t  t  p 

k 

Fig.  A1 ,  4  .01  r—  4.01  radians  angular  width,  with  magnitude  bounds 
2_ 

k  k  ' 


.99  A  ,  1.01  A 
.  k2  kl 


.  Let  A  in  Fig.  A1  be  a  trial  value  of  £  ,  so 

o 


P  P 

u  rO 

—  +  £q  is  the  vector  originating  at  point  —  in  V  and  terminating 

at  A.  Bounds  on  £  may  be  obtained  so  that 
0 


sup 

p 

P  0 

inf 

P„ 

P 

—  +  £ 

P  0 

satisfies  A2(/,2)and  A2(3)  in  It=  “). 


It  is  easily  seen  that  a  very  conservative  boundary  for  £q  in 
Ij.  is  the  vertical  line  s  =  -o,  with 


o  =  a 


_0  [-99A  -  l.oi 
1  "  2  I  Xk2  -  k1 


>  0 


i.e.,  £q  on  the  right  of  the  line  s  =  -a,  satisfies  A2(l,  2)  in  1^. 


(A3) 


(A4a) 


p.  1  _ 

ro 

1  +  pq  1 

p 

—  +  £ 

P  o 

£  B1 1  p  i ,  or 


Ho/  l‘ 
P  o 


This  is  easily  satis¬ 


fied  if  the  above 


(A4b) 


2  B 


1 


Therefore,  choose 

a  =  smaller  of  (o^,  a  )  (Me) 

Thus,  the  problem  is  to  find  £o(s)  such  that  |fco(jw)|  is  outside  the 
circle  C  in  Fig.  A1  for  and  to  the  right  of  the  line  s  ~-a  for 
w  >  It  is  obviously  very  easy  to  find  such  an  £Q(s)  which  also 

has  all  its  poles  and  zeros  in  the  interior  of  the  left  half-plane, 
with  any  desired  finite  excess  of  poles  over  zeros,  and  which  further¬ 
more  has  the  property  shown  in  Fig.  Al,  i.e.,  lies  on  the  right  of 
Ss-o,  for  all  a).  For  example,  let 


2  i 


A-  = 


ox 


(i  +  -5-)  n  0  +  } 

V  i  ui 


e  any  desired  finite  number,  = 


21  iii. 

larger  of  (10u>t,  — — -),  w-+1  =  100w. 


Note  that  it  would  be  impossible  to  guarantee  the  existence  of  the  desired 
iQ  if  pQ  was  nonminimum-phase  (Horowitz  and  Sidi  1978). 
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It  is  conceivable  that  even  though  A2(l)-A2(3)  are  satisfied, 

p 

A2(4)  is  not  satisfied.  Consider  the  zeros  of  1  +  pg  -  1  +  z  —  or 

o  p 

_  *p 

of  (“- +  «.Q).  Recall  that  in  [0,  lies  in  the  circle  of  radius 

Y  which  is  inside  the  circle  C  of  Fig.  Al,  while  Z  lies  outside  the 
L  o 

larger  circle  C.  In  °°),  *  iies  on  the  right  of  the  line  s=-  o 

while  is  contained  in  V  of  Fig.  Al.  Also,  ^(jw)  lies  on  the 

right  of  s  =  -a,  Vu>.  Hence,  the  vector  ^  +  z  does  not  encircle  the 

p  o  - 

origin  clockwise  (or  alternatively  does  not  so  encircle  VpcP, 
and  the  system  is  stable. 

Application  to  mio  system 

In  the  mio  system  (4a),  the  loop  function  *U(j  =  9UU/PUU  must  handle 

the  n  tuv  problems  v  =  l,2,...,n.  The  bounds  on  S,uuQ  will  be,  in  general, 

different  for  each  v  with  its  own  toxv,  Itv  of  (A3)  and  oy  of  (A4c). 

Let  Iu  =  max  Itv,  £QU  =  max  £QXv,  oy  =  min  oy  be  the  design  parameters 

for  £UU().  Obviously  such  a  S,uuo  is  satisfactory  for  all  n  t  problems. 

In  Sec  3  (just  before  3.1),  there  was  noted  the  desirability  of 

satisfying  (6)  on  the  boundary  of  a  semicircle  of  radius  in  the  right 

half-plane.  This  requires,  in  addition  to  the  previous,  rewriting 

i6 

Sec  A3,  replacing  ju>  by  u>HeJ  ,  6e  (0,  ti/2].  The  development  is  easier  if 

is  large  enough  so  that  each  p.j  =  k^/s  1J  on  s  =  w^e'19  .  Clearly, 

there  will  emerge  bounds  on  on  s  =  wueJ  ,  which  will  have  to  be 

uuu  n 

satisfied,  in  addition  to  those  on  s  =  jw.  Obviously,  such  bounds  can 
always  be  satisfied  by  suitable  shaping  of  *uu0  ,  so  that  I Auuo I  1s 
large  enough  on  • 


There  remains  the  design  of  f(s),  inasmuch  as  ^(s)  =  g(s)  pQ(s) 
only  determines  g(s).  Mote  that  only  guarantees  that 


sup  lPo 


?  lp  +£°l  <BMof 


inf 

P 


£♦1 

P  o 


A(oj) 


A2(l); 


f(s)  is  chosen  so  that  |t(joj)|e  {A(u>),  B (u>) ] .  For  example,  suppose 
A(u>-|  )  =  .9,  B (to^ )  =  1.04,  and  at  u^,  Sjv^|~p^  +  =  100,  while  inf 


V* 

P  ° 
f*. 


=  90  with  K0(jio^)|  =  80.  The  range  of  |t(ju^)j  = 


is  therefore  {.8| f (jt^ ) | ,  .889 j  f ( jcoj ) |] ,  so  we  need 

|>.9,  .889|f(ju^ ) J <1 . 04 ,  giving  the  permissible  range  of 
for  |f(jwj)|.  In  this  way,  the  bounds  on  |f(jw)|  are 
found  and  it  is  always  possible  to  find  an  f (s )  with  left  half-plane 
poles  and  zeros  which  satisfies  such  bounds. 

The  above  procedure  in  all  its  details,  is  not  recommended  as  a 
practical  design  procedure.  Simplifications  were  made  to  make  the 
proof  easier,  but  the  loop  bandwidth  is  much  larger  than  necessary. 
Its  primary  purpose  is  as  an  existence  theorem.  A  practical  optimum 
design  procedure  based  on  these  ideas,  but  without  the  rigor,  has 
been  given  in  (Horowitz  and  Si di  1972)  and  used  a  great  deal  with 


—  +  l 
P  o 


•8|f(jwj) 
9  1.04 


8’  .889 


considerable  success. 
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A5.  Extensions 

(1)  It  is  possible  to  have  A(w)  =  8(10),  Q(w)  =  0  in  A2 ( 1 ) ,  (3) 
at  a  finite  number  of  u  values,  by  choosing  g  infinite  at  these  points. 
The  sensitivity  zeros  can  be  single  or  multiple. 

(2)  Some  or  all  pePcan  have  zerox  on  the  jw  axis.  If  these  zeros 
are  precisely  known  (unlikely),  g(s)  can  be  assigned  poles  there.  Other¬ 
wise,  t(ju)  must  be  zero  at  these  points  for  such  p,  requiring  obviously 
much  more  careful  statement  of  the  tolerance  on  t(ju)  and  trf ( joi)  near 
such  points. 

(3)  The  most  significant  extension  is  that  Constraint  Al(l)  can 
be  relaxed.  There  can  be  uncertainty  in  the  order  of  the  plant  due  to 
disappearing  poles  and  zeros--closely  related  to  the  problem  of  singu¬ 
lar  perturbations  (Porter  and  Tsingas  1978). 


A6.  Disappearing  Poles  and  Zeros 
Let 


P 


m 

p^l  (1+sa . ) 
n  (l+sbj) 


=  P1  'P(s) 


(A5) 


with  a.elO,  a.  1,  b.e[0,  b.  Jand  p.e P satisfying  A1 .  The  question  of 

1  IX  J  JX  1 

concern  is:  "For  what  m,  n  values  can  the  loop  transmission  be  arbi¬ 
trarily  large  over  an  arbitrarily  large  bandwidth  but  still  be  practical, 
i.e. ,  go  to  zero  as  w  *■  «  with  any  desired  finite  excess  of  poles  over 
zeros?"  For  such  m,  n  any  tolerances  satisfying  A2  but  otherwise  ar¬ 
bitrary,  can  be  satisfied. 


m 


In  Fig.  2  and  using  (A5) 


t(M  =  las _ =  fSPll _ 

K  ’  1  +  gp  1  +  gp,v 


ft 

1  +  z  ' 


(A6) 


The  question  posed  can  be  answered  by  referring  to  the  logarithmic 
complex  plane  (Nichols  chart)  in  Fig.  A2.  The  intersections  of  the 
zero  db  line  with  the  vertical  lines  (2n+l)  180°,  n  =  0  ±  1,  .  .  . 
is  the  point  -1.  Because  of  uncertainty,  Z  =  gp^V  is  not  a  point  but 
a  set  {£}  denoted  here  as  the  template  of  Z,  ^(w)  which  occupies 
some  region  in  the  complex  plane--Fig.  A2.  The  shape  of  is 

•X 

that  of  {pj  because  there  is  no  uncertainty  in  g.  The  latter 
permits  the  translation  (but  not  rotation  of  £)  in  the  complex 
logarithmic  plane,  horizontally  by  arg  g  and  vertically  by  |g|  in 
db.  For  some  finite  w  range,  large  jfcj  is  needed  so  lies  above 
the  zero  db  line,  e.g.,  (co^ )  in  Fig.  A3.  At  large  enough  u>,  |t| 
must  be  very  small  (-►O  asur*0)  so  must  be  well  below  the  zero  db 
lines  and  continue  downwards  to  -  ».  In  the  transition  of  ^from 
above  to  below  the  zero  db  line,  it  must  not  interest  -1,  nor  en¬ 
circle  it.  Hence,  the  width  of  must  be  restricted  to  <360° 
for  some  u>  interval  in  which  £^can  squeeze  in  between  two  -1  points 
on  Its  way  downward  (Fig.  A£).  But  we  want  arbitrary  sensitivity 
reduction  for  arbitrary  bandwidth.  This  requires  3  5  f°r 


m  >  uj|_|  the  width  of  ^j(u)  <  360°. 


37 


Consider  now  the  shape  of  ^  which  is  that  of^  which  is 
that  of £7p1  4^  in  the  Nichols  chart.  Constraint  A1  (3)  assures 
that  at  large  enough  w  the  width  of  -*■  zero  degrees,  so  it 

is  entirely  a  question  of  the  width  of  .  Consider  any  factor 


of  Y,  e.g. , 


1 


1  +  job 


,  be  [0,  m].  At  any  specific  o,  the 


maximum 


width  of  the  template  is  tan"1  wb^  =  tan”1  am.  For  tun  =  20,  the 
template  is  OAB  in  Fig.  A3  (O  at  b  «*  0,  B  at  ob  =  20);  for  urn  =  100 
it  is  0ABU  and  for  urn  =  1000,  it  is  0ABUV.  For  two  such  independent 
factors,  it  is  easy  to  find  the  new  template.  This  is  done  in  Fig.  A3 
for  =  20.  Draw  0"A"B"  =  OAB  and  position  0“  at  points  on 

OAB  (because  of  the  independent  uncertainties).  The  result  is 
ABEJDEFO  = 

”) ^T+jojb  ) ^  *  ^1’  ^2  temP^ate  a 

zero  factor  (1+jwa)  is  obtained  by  reflecting  that  of  the  pole  if 

(“^max  =  (wb)  max*  Saving  OA'B'  in  Fig.  3. 

As  w  increases,  the  contribution  of  each  factor  (pole  or  zero)  -+•  90° 
in  width.  Therefore,  while  theoretically  four  factors  can  be  admitted 
between  two -1  points,  stability  margins  dictate  a  maximum  of  three. 

But  this  is  only  a  necessary  condition,  because  it  must  also  be  possible 
to  decrease  |L|  from  arbitrary  large  to  small  values.  This  basically 
means  that  over  arbitrary  large  frequency  range  Arg  L  must  be  over¬ 
whelmingly  negative.  The  extreme  right  side  of ^  must  then  lie  on  the 
left  of  the  0°  line  in  Fig.  A2.  One  disappearing  pole  or  zero  poses  no 

problem,  but  two  do  because  the  left  side  of  J  will  then  intersect  -1 

Jt 


be  made  rigorous,  of  course.)  Note  that  this  is  for  the  most  demanding 
situation.  If  only  stability  is  desired,  then  any  finite  number  of 
disappearing  poles  or  zeros  can  be  handled. 


APPENDIX  2 

A7.  CONDITIONS  THAT  $  MAPS  H(n2)  INTO  ITSELF 

The  constraints  on  Huv  were  given  in  2.1.  The  first  is  satisfied 
if  (6)  is.  From  Appendix  1,  it  is  seen  that  at  finite  to,  the  only 
possible  difficulty  is,  if  in  Fig.  2,  duv ( Juj)  is  unbounded  but  9ul|Puve 
is  bounded.  If  such  unboundedness  of  duy  is  at  a  finite  number  of  u> 
values,  it  is  possible  to  assign  poles  there  to  gu(J.  Since  dyv  is  a 
function  of  the  plant  parameters  (Eq.  4b),  an  infinite  number  of  such 
poles  is  conceivable,  in  fact  is  so  in  practice  if  there  is  one, 
because  of  inevitable  uncertainties.  How  can  duy  be  unboundeo  while 
Puve  is  bounded?  From  (3,5) 


i?u 


1U  IV 


A  -  Y  (-1 )  p .  A .  +  p  A 

/•  '  '  Hiu  ni  Huu  uu 

ij*u 


(A7) 


Hence,  such  a  situation  is  possible,  if  A^u  has  such  a  pole  which  is 
cancelled  by  a  zero  of  p^  and  is  not  present  in  PUUAUU  e.g.  n  =  2, 

P2j  has  a  pole  at  ±ja,  p^  has  a  zero  there  and  P^Pgg  does  not  have 
a  pole  there.  Such  situations  are  therefore  not  allowed.  It  can  be 
argued  that  they  are  in  practice  impossible  (Horowitz  and  Shaked  1975), 
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because  P]2>  P21  must  involve  different  physical  parameters  and  the 
uncertainties  in  each  cannot  be  100%  correlated.  Irrespective  of  thi 
argument,  the  constraint  is 


I ^jj(>3u))|  is  uni  fonnly  bounded  on  [O,10)  and  M,  V  i^j 


(A8) 


The  problem  is  more  difficult  as  w  •>  »  and  is  treated  in  Section  3. 
To  satisfy  constraint  (2)  on  in  2.1,  note  first  that  if 
:(jw)  =  Zj  +  jz2  with  z]t  z2  real  functions  of  w,  then  Jz|*  is  bounded 


if  z'  is  bounded.  For 


-  <  n  *  -lilLLv 


1  1 


2  2 
Z1  +  Z2 


1 

—  If  z1  =  z|  +  jz\ 


is  bounded,  so  are  z^,  z ^  and  \z\  1 .  Consider  separately  t1  (joj)  and 
^Tduv  duv^  of  The  former  can  be  written 


t  1  =  (JL.V  =  f  +  fl32L±3lJ>± 

"l  (1*9P)2 


1  =  9P-  O’  tw  S  *  9UU 


(A9) 


of  Fig.  2.  Since  f£/(l  +  £)  is  bounded  by  the  appropriate  Buv>  the 
first  term  on  the  right  of  (A9)  needs  only  uniform  boundedness  o*  f ' , 
which  is  easy  as  f  is  chosen  by  the  designer.  Obviously,  the  only 
possible  difficulties  with  the  second  term  of  (A9)  are  g’,  p'  which 
can  be  infinite  only  at  jw  axis  poles.  However,  at  such  poles,  the 
denominator  forces  the  second  term  to  be  bounded.  At  large  w  where 
gp,  g',  p'  ->  0,  there  is  obviously  no  problem. 


Next,  consider  (Ttjuvduv^0^  writtfn1  as 


(A10) 


y‘ 


A'.s. .  + 

1  +  9P  (l  +  gp)2 


g,p  of  (A9) ,  d  -  duv  . 

From  (A8),  d  and,  sj ruse. 4  j,S.  3  rational  function,  d',,  arg  jjnit;',mly 
bounded  on  the  jto  axis.  From  the  previous  discussion  re  p',  g1  etc.  and 
Appendix  1,  y'  can  be  uniformly  bounded  on  any  finite  o>  range  by  proper 
choice  of  g.  At  large  w  each  of  d,  d',  p,  p*  etc  0. 

It  is  easy  to  see  that  constraint  (3)  is  satisfied,  because  it  has 
been  required  that  the.  elements  of  P,  G,  F  all  -*•  0  as  s  -*■  °°. 
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Fig.  1  Multiple  input-output  two  matrix  degree-of-freedom  feedback 
structure  c  =  Tr,  T  -  [t- -].  c  =  [c,  ...  cj\  r  =  [r,  ... 
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Fig.  5.  "Equivalent”  frequency  domain  bounds  and 
experimental  results  for  various  plant 
parameter  sets 


